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ABSTRACT 

We report the results of a cosmic shear survey using the 4.2m William Herschel Telescope on 
La Palma, to a depth of R = 25.8 {z rs 0.8), over4 square degrees. The shear correlation func- 
tions are measured on scales from 1' to 15', and are used to constrain cosmological parameters. 
We ensure that our measurements are free from instrumental systematics by performing a se- 
ries of tests, including a decomposition of the signal into E- and _B-modes. We also reanalyse 
the data independently, using the shear measurement pipeline developed for the COMBO- 17 
survey. This confirms our results and also highlights various effects introduced by different 
implementations of the basic "KSB" shear measurement method. We find that the normalisa- 
tion of the matter power spectrum on 8 ft, _1 Mpc scales is er 8 = (1.02 ± 0.15)(0.3/O m ) ' 5 , 
where the 68%CL error includes noise, sample variance, covariance between angular scales, 
systematic effects, redshift uncertainty and marginalisation over other parameters. We com- 
pare these results with other cosmic shear surveys and with recent constraints from the WMAP 
experiment. 
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1 INTRODUCTION 

Weak gravitational lensing by large-scale structure, or "cos- 
mic shear", has emerged as a powerful cosmological probe, 
as it is directly sensitive to foreground mass (f o r review s, see 



Bartelmann & Schneiderl l200fi lBernardeaul ll999t iMellied fl999l 
02) 



Refregier 2003; Wittman 2002). A measurement of cosmic shear 
is therefore closely tied to cosmological theories, which are prin- 
cipally concerned with the distribution of dark matter. In particu- 
lar, the systematic biases of this technique are not limited by un- 
known physics such as biasing l Dekel & Lahav 1999; Grav et al. 
[ 20021: iHoekstra et all l2002bt ISmith et alJ 120031: IWeinberg et alj 
2000) or th e mass-temperature relation for X-ray selected galaxy 
clusters iHuterer & White! 120031; IPierpaoli. Scott & White! 1200 it 
IViana. Nichol & Liddlel2002l) . 

Cosmic shea r surveys ar e rap idl y growing in size 

and precisio n (Bacon et al.l | 2003l: iBrown et alj 12003 

Haman^eUilJj20 i 03^_jHoe IstraeTal] 
Refregier Rhodes. & Groth 



2002; 



2002a i: Ijarvis et alj 12 003 
Rhodes, Refregier & Groth 



200ltlvan Waerbeke et alj I2OO2I) . Cosmological parameter con- 
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straints from these surveys are now approaching the precision of 
other methods. 

However, cosmic shear surveys can be subject to several sys- 
tematic biases of their own. Imperfect telescope tracking, telescope 
flexure or optical misalignment within the camera, even at a level 
that is acceptable for most purposes, can artificially distort images 
and align the shapes of distant galaxies in a way that mimics cosmic 
shear. 

The survey described in this paper represents a culmination of 
effort at the William Herschel Telescope. We have combined the 
experience of instrumentalists with detailed image simulations and 
careful data analysis to control th e various sources of syst ematic 
error. Our first cosmic shear paper I Ba con. Refregier & Ellisl2c5ocl) 
reported an initial detection of cosmic shear using a 0.5 square de- 
gree surve y with the William Herschel Telescope (WHT). The sec- 
ond paper iBacon et all2003l) compared the WHT shear signal with 
an independent measurement using the Keck II telescope, and ex- 
amined systematics from these two very different instruments. In 
this paper, we extend our WHT survey to cover 4 square degrees 
to constrain cosmological parameters, while paying great care in 
monitoring and correcting systematic effects. We also now test our 
entire pipeline for shear measurement and cosmological parameter 
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estimation by comparing it to extern al code, developed indepen- 
dently for the COMBO- 17 survey bv lBrown et"ail <2003l) . 

This paper is organised as follows. In fj2]we describe our sur- 
vey strategy and observational parameters. In J3|we present our 
results and draw constraints upon cosmological parameters. In i]4] 
we test for the absence of any systematic errors. In f|5] we present a 
second set of results obtined via an independent shear measurement 
pipeline. We conclude in J6| 



2 OBSERVATIONS AND DATA REDUCTION 

We have acquired 4 square degrees of imaging to R — 25.8 (for a 
point source at 5a) with the Prime Focus Imaging Camera (PFIC) 
of the William Herschel Telescope on La Palma. The median see- 
ing was 0.69" and no exposures had seeing worse than 1". The 
pixel size is 0.24". As shown in figureQ pointings were scattered 
randomly in a pencil-beam survey between galactic latitudes of 30° 
and 70°. This was tuned to provide ~ 1.5 stars per arcmin 2 , with 
which we could measure the Point Spread Function (PSF) across 
each field. The only selection criterion was to avoid foreground 
stars brighter than R ~ 11 in the Digitised Sky Survey or APM 
(Automated Plate Measuring machine) catalogues. 

Cosmic shear statistics have already been presented from the 
first square degree of this survey in lBacon e t al. (2003). That data 
consisted of eight 8' x 16' images and eleven 16' x 16' images taken 
after the addition of a second, identical CCD to the PFIC. During 
June and August 2002, we observed an additional 41 16' x 16' 
pointings. When combining statistical measurements from all of 
these fields, we weight the contribution of the large fields to be 
twice that of the small fields. This scheme is ideal for measure- 
ments on small angular scales, where the large fields do indeed 
contain twice as many pairs of galaxies as the small fields. The 
scheme is non-optimal in exploiting the extra signal on large scales 
contained via additional pairs across different CCDs. However, in 
general, as both large and small fields were randomly deployed, no 
bias should result from our weighting scheme. 

For each survey field we took four 900s exposures, each 
dithered by a few arcseconds from the last. This strategy enabled 
a continual monitoring of astrometric distortions within the tele- 
scope, cosmic ray removal, and lower overheads in the event of in- 
clement weather. Dat a reduction then proceeded for the exposures 
exactly as in lBacon et alj ^003): see that paper for more details. 
After bias subtraction and flat fielding, fringing remained in the 
i?-band images. This could have been prevented by observing at a 
shorter wavelength, but at a cost to the observed number density 
of background sources. To remove this, a fringe frame was com- 
piled from all the exposures in each night. A multiple of this was 
subtracted from each image which minimised fringing, to a neg- 
ligible level < 0.05% of the background noise. The four dithers 
for each field were then re-aligned (using linear interpolation be- 
tween adjacent pixels to allow sub-pixel offsets) and stacked (with 
3cr-clipping to remove cosmic rays). 

Objects were located on the final images using HHNDPEAKS 
from the IMC AT package bv lKaiser. Sau res & Broadhurst| i 1 9951 
KSB). Following the recommendations of M^ssev^t^m200ll) . ob- 
jects within 10" of saturated stars or 5" of the edge of the CCDs 
were masked and removed. We also removed noisy objects from 
the catalogue with cuts in size, signal-to-noise and ellipticity of 
r 9 > 1 pixel, v > 1 5 and |e| < 0.5. This is the same proce- 
dure to that adopted in iBacon et alj l2003h . and leaves 15.2 galax- 
ies per arcmin 2 in the final catalogue, with a median magnitude of 
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Figure 1. The locations of WHT pointings in our deep pencil-beam survey 
strategy. The galactic latitudes were tuned to provide sufficient foreground 
stars within each image to successfully model and correct for variations in 
the PSF. 



23.5 ± 0.2. According to lCohen etai] fcOOOl) as before, this corre- 
sponds to a median source redshift of z s ~ 0.80 ± 0.06. 

Galaxies' observed ellipticities, e, were formed from combi- 
nations of their Gaussian-weighted quadrupole moments, then cor- 
rected for convolution with the telescope's PSF. The shape of the 
PSF was measured from stars in each image, then interpolated to 
the positions of galaxies via a polynomial fit within each CCD sep- 
arately. The order of this polynomial was varied from field to field: 
it was occasionally raised to fourth order but lowered to second or- 
der where possible. This limited spurious power on small scales and 
particularly around the edges of the field. Even so, an acceptable 
polynomial fit to the PSF was not always possible. We conserva- 
tively discarded 14 fields, all of which were otherwise within spec- 
ifications, but whose stellar ellipticities remained correlated after 
correction. More sophisticated PSF fitting met hods have re cently 
been suggested by Hoekstr J 120041) and ljarvis & Jainl i2005;). Such 
techniques may improve the interpolation, and thus increase the 
number of useable fields, but they have not been investigated here. 

The shear susceptibility factor for each galaxy, P 1 , was de- 
termined from the higher-order shape moments of the galaxy and 
the interpolated PSF. The shear susceptibility is a notoriously noisy 
quantity; in this implementation of KSB, we chose to fit P 7 with 
a cubic polynomial as a function of galaxy size, r B . We can then 
form shear estimators for each galaxy, 7 = e/P~ r . We apply a final 
calibration factor of (0.85 ± 0.04) ~ to these she ars. This calibra- 
tion factor was determined by B acon et alj feOOll) . using our shear 
measurement pipeline upon simulated WHT images with a known 
input signal. 

Each set of four dithered exposures were also used to con- 
tinually mo nitor astrometric d istortions within the telescope. As 
observed in Bacon et al. (2001), these closely follow the engineer- 
ing predictions in the PFIC manual of 7tan g cntiai = 0, 7 ra dial = 
—8.2 x 10 _5 r 2 with r measured in arcminutes from the field cen- 
tre. This affects the shear estimate of an average galaxy in a large 
exposure by less than 2 x 10~ 3 . A galaxy in the corner of a large 
exposure is altered by 1 x 10~ 2 . Since this is no longer an entirely 
negligible effect in our enlarged survey, we subtract this distortion 
from the final shear catalogues using the shear addition and sub- 
traction operators in Berns tein & Jarvisll2002l) . 

The resulting distribution of shear estimators is shown in fig- 
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Figure 2. The distribution of measured shear estimators from all galax- 
ies used in the WHT survey. A fitting function is given in the text. The 
two components of shear are shown superimposed; their distibutions are 
indistinguishable. Notice, however, the extended wings typical of the KSB 
method. 



ure|2] Both components of shear are well-fit by the normalised 
proabability density function 



P( 7i ) w 0.26 



exp(-:r 2 /0.65 2 ) 
(x 2 + 0.27 2 ) 



(1) 



Extended, Cauchy-like wings are a salient feature of the KSB 
method, where ellipticity measurements are formed from ratios of 
(noisy) quadrupole moments. 



3 RESULTS 

3.1 Shear-shear correlation functions 

The power spectrum of the weak lensing shear is given by 



C7 = 



_9_ (Ho 
16 V c 



<>- 



g(x) 



Da{ X ) 



(-,*) dx, 



(2) 



where \ is a comoving distance; \h is the horizon distance; Da (x) 
is an angular diameter distance; g(x) is the lensing weight function; 
and P(k, z) is the underlying 3D distribution of mass in the uni- 
verse. The two-point shear correlations functions can be expressed 
in terms of the power spectrum as 
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1 

4tt 
1 

4tt 
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(3) 
(4) 



These can be measured by averaging over galaxy pairs, as 

cm = (Yii^Yiir + e)) (5) 
02(0) = (72M 72 r (r + e)), (6) 

where 9 is the separation between the galaxies and the superscript r 
denotes components of shear rotated so that 71 (72) in the first 
galaxy points along (at 45° from) the vector between the pair. A 
third shear-shear correlation function can be formed. 



C^) = (7i»72> + 0)> + <72(r)7i r (r-+0)) 



(7) 



for which the parity invariance of the universe requires a zero sig- 
nal. Cz(ff) can therefore be used as a first test for the presence of 
systematic errors in our measurement. 






Figure 3. Correlation functions of the shear field measured in our 4 square 
degree WHT survey. The solid data points show our measurement. The in- 
ner error bars are for statistical errors only; the outer error bars also include 
full non-Gaussian sample variance. The blue circles show measurements 
of our data by the COMBO-17 pipeline, artificially adjusted to correct for 
its slightly higher source redshift distribution and permit easy comparison. 
The solid line shows the theoretical prediction for a ACDM model with 
f^m = 0.3, Q\=0J, r = 0.21 and as = 1.0, assuming a median source 
redshift of z s = 0.8 and using the fitting functions of Smit h* et alj i20o3) . 
The dotted lines show similar theoretical predictions, but with eg rang- 
ing from 0.7 (bottom) to 1.2 (top). The dashed lines barely visible above 
the x axis show the correlation of galaxy shears with the observed PSF 
anisotropy. 



To perform the measurement in practice, we first measure 
the shear correlation functions in each field, using all galaxy pairs 
within various 9 bins. To obtain a combined result for the entire sur- 
vey, we then average the binned values for each field. We find that 
the measured correlation functions are quite sensitive to changes 
in the binning scheme. The highly non-Gaussian distribution of 
shear estimators (see figure[2j, plus additional non-Gaussianity in- 
troduced by computing correlation functions, means that the in- 
clusion in one bin of a small number of outlying fields can sig- 
nificantly bias measurements. Both the correlation functions and 
the subsequent constraints on cosmological parameter constraints 
move within their full la statistical error bars. During the aver- 
aging, we therefore introduce 3u-clipping in each bin to remove 
some outliers, and further 3<r-clipping in C3 (9) and the star-galaxy 
cross-correlation functions Cf G (9) + Cf G (#) to eliminate flaws 
with PSF correction (see £14.21 . This leaves between Nf — 40 and 
Nf = 43 fields used for each angular bin. Our final choice of bin 
size yields representative central values. The result is shown in fig- 
ure[5j We shall separately estimate the amount of instability caused 
by binning, and increase the error upon cosmological parameter 
constraints accordingly. 

In order to derive any constraints on cosmological parameters, 
it will also be necessary to know the covariance of d (9) between 
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Figure 4. Covariance matrix of different angular bins of the shear cor- 
relation functions shown in figure |3] from small to large 9, then from 
small to large 9 again. For example, the bottom-left quarter shows 
cov[Ci(0),Ci(#)]. The diagonal elements within this are the variances 
in each bin. 



different angular bins. Our pencil-beam survey strategy with many 
independent fields makes it easy to measure their covariance matrix 



cov[C , i (0),C' j (0)] 



(8) 



/=i 



where the summation is over all fields, and the superscript ' de- 
notes correlation functions calculated in one field alone. This ma- 
trix is depicted in figure [4] and shows the significant covariance, 
especially between adjacent bins. 



3.2 Cosmological parameter constraints 

We now use a Maximum Likelihood method to determine the con- 
straints set by our observations upon the cosmological parameters 
il m , the total mass-density of the universe, and as, the normali- 
sation of the matter power spectrum at 8 h^ 1 Mpc. The analysis 
directly uses the obse rved correlation functions C\(6) and dip), 
proceeding as in IB aeon et alJ 120031) . except that theoretical pre- 
dictions for the non-linear powe r spectrum a re calculated via the 
updated fitting fun ctions of lSmifh et alJ 120031) rather than those by 
iPeacock & Doddsl ll996l) . This has the effect of lowering our fi- 
nal constraint on <7gf2^ by about 5%. We use the fitting functions 
for the linear transfer function suggested by iBardeen et all 1983) . 
Note that, although we shall perform an E/B decomposition, we 
fit Ci and C2 rather than the E mode signal to avoid degeneracies 
arising from our survey's finite size. We will use the E/B decom- 
position as an a posteriori consistency check for systematics on rel- 
evant scales; because of such contamination, we discarded the first 
and last data points shown in figure|3|(see discussion in *|4j. 



1.5 



1.0 



0.5 




0.2 



0.4 



0.6 



0.8 



1.0 



Figure 5. Constraints upon cosmological parameters Q m and erg, from a 
maximum-likelihood analysis of our WHT cosmic shear survey data. The 
68.3% (solid), 95.4% (dashed) and 99.7% (dotted) confidence limits include 
statistical errors and non-Gaussian cosmic variance. However, they include 
neither the calibration of the shear measurement method, nor uncertainty in 
the source galaxy redshift distribution. These sources of error are consid- 
ered separately in the text. 



The theoretical correlation functions were first calculated 
from equation on a 2D grid across the 0, m vs as plane. The 
power spectrum shape parameter was set to T = 0.21, consistent 
with re cent observat ions of clust ering in galaxy redshift surveys 
IPercival et al. 2001; Szalav et al. 2003). The median redshift for 
source galaxies was fixed to z s = 0.8 for WHT and z s — 1.0 for 
Keck. Errors on these parameters will be propagated separately into 
our final constraints. 

We then fitted the observed shear correlation functions d{9) 
to the theoretical predictions calculated at the centres of each bin 
t(6), computing the log-likelihood function 

v t 



X = 



tie., Q r . 



*>)) 

cov[Ci(9),Cj(#)] 



(d(ti)-t(d,£l m ,(T S j) (9) 



throughout the grid. We thus explore parameter space in this plane, 
and minimise \ 2 to find the best-fit cosmological model. To com- 
pute confidence contours, we numerically integrate the likelihood 
function 



,0-8) 



-X 2 /2 



(10) 



Note that the matrix inversion in equation {5J is unstable, 
possibly for the same reasons as the binning instability. Since 
cov[Ci(#), C2{0)] is constitent with zero but noisy, we solve this 
problem in practice by forcing it to zero. The inversion is then well- 
behaved. 

Our constraints on cosmlogical parameters are presented in 
figure |5| and the constraints from our Keck survey iBacon et alJ 
2003) are reproduced in the same format in figure|6| In both cases, 
the contours show 68.3%, 95.4% and 99.7% confidence limits, in- 
cluding statistical errors and non-Gaussian sample variance. They 
reveal the well-known degeneracy between Q m and as when using 
only two-point statistics. 
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Figure 6. Constraints u pon cosmol ogical parameters from the Keck cosmic 
shear survey bv lBacon eta?] 1200 31 showing the 68.3%, 95.4% and 99.7% 
confidence limits as in figure |5] Only one edge of the 99.7% confidence 
contour is visible inside this parameter range. 



A good fit to the 68.3% confidence level from our WHT data 
is given by 

, O \ 0.50 



ni m y 

Vo.3 7 



1.02 ±0.13, 



while the Keck data is well-fit by 

, n \ 0.52 
1 £t? 



08 



/ i 2 m 
Va3 



= 1.01 ±0.19, 



(11) 



(12) 



The multiplication of the respective likelihood functions provides 
a constraint from a combined survey. Such confidence contours are 
shown in figure^ with the 68.3% confidence level well-fit by 



I 0.3/ 



= 1.02 ±0.12 , 



(13) 



for 0.1 < fi m < 0.7. 

Note that all of these constraints include only the statistical 
error and sample variance. We can propagate other sources of error 
by noting that 

r< Cc'\ o!-46 2.45 L65ti-0.11 / D 7\-2 t ,.^ 

Ci(5)ocS2 m er g z s F (P'J , (14) 

where i = 1 and 2 and P 1 is the shear calibration factor, in a 
fiducial ACDM cosmological model with Q, m = 0.3, f^A=0.7, F — 
0.21 and erg = 1.0. Adding in turn to our constraint <13> : a 6% 
additional binning instability, a 10% source redshift uncertainty, a 
15% prior on F, and a 5% shear calibration uncertainty, gives a final 
68.3% CL constraint for the combined survey of 

'ft ^°' 5 



/ftnV 
V 0.3 ^ 



= 1.02 ± 0.12 ± 0.06 ± 0.066 ± 0.006 ± 0.04 (15) 
1.09 ±0.15, (16) 



where the various errors have been combined in quadrature on the 
second line. This result now includes all contributions to the total 
error budget: (non-Gaussian) statistical noise and sample variance, 
covariance between different angular scales, binning instability, 
shear calibration error, source redshift uncertainty, and marginal- 
isation over F. 
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Figure 7. Constraints upon cosmological parameters for the combination 
of both WHT and Keck surveys, showing the 68.3%, 95.4% and 99.7% 
confidence limits, as in figurel5l 



3.3 Shear variance 

For historical reasons, cosmic shear results are often expressed as 
the variance of the shear field in circular cells on the sky. For a top- 
hat cell of radius 9, this measure is related to the shear correlation 
functions by 



< Iff > 



CJ{1) [^(tO)] 2 IAI (17) 



(18) 



where we have used a small angle approximation involving the 
Bessel functions. Note that the shear variance is more strongly cor- 
related on different angular scales in this form than they are as cor- 
relation functions. 

In practice, data is not available on all the scales necessary to 
perform this integration. The correlation functions have not been 
calculated on scales smaller than 1', and are contaminated by sys- 
tematics on scales smaller than 2' (see j4j. We determine the deficit 
in the measured values as a function of 8 by extrapolating the 
data through these scales using the theoretical predictions given 
by the best-fit cosmological model determined in A3. 21 This deficit 
(~2x 10~ 5 at 3' and 1 x 10~ s at 5') is then added back on to our 
measured data points. 

We present our results as the variance of shear in cells, 
and compare them to those from similarly deep lensing sur- 
veys in figure [8] These surveys use data from 8.5 square de- 
gree, z 3 = 0.8 VIRMOS-Descart survey on the 3.6m CFHT by 
IVan Waerbeke. Mellier & Hoekstrd {2005 ) ; the 1.25 square degree, 



z s = 0.85 COMBO -17 survey on the 2.2m La Silla telescope by 
iBrown et alj 1200 3l) : the 0.36 square degree, z s = 0.9 Medium 
Deep Survey (MDS) with the Wide F ield an d Planetary Camera on 
HSTby Refreaier Rhodes, & Groth (2002); 0.27 square degrees of 
random fields observed to z s = 1 .0 in parallel-mode with the Space 
Telesc ope Imaging Spectrograph (STIS) on HST by Rhod eset all 
(2004); and the 0.6 square degree, z s = 1.0 pencil-be am survey 
using the 10m Keck II telescope bv lBacon et all 120031) . For ease 



6 R. Massey et al. 



EST MDS Refregier et al. (2002) : 
■ Keck, Bacon et. al. (2003) 1 
X COMBO- 17, Brown ot al. (2003) 
A HST STIS parallels, Rhodes et al. (2004) : 
□ VIRMOS-Deseart, Van Waerbeke et al. (2004) : 
♦ WHT, this work 1 




Figure 8. Shear variance in (circular, top-hat) cells, as a function of the 
radius of the cells. Our results are compared to those from similarly deep 
surveys by other groups. For ease of comparison, all results have been 
rescaled as if their median source redshift had been the same (see text). As 
in figure[3] the theoretical curves are predictions for a ACDM cosmology, 
with Q. m = 0.3, S1a = 0.7, T = 0.21 and <rg ranging from 0.7 (bottom) 
to 1.2 (top). 



of comparison between these different surveys, all the results have 
been scaled to the values that would have been obtained if their me- 
dian source redshift had been z s — 0.8, and using the theoretical 
prediction that oc zl' 65 . Results from deeper surveys are thus 
shown slightly lower here than in their original papers. 



4 TESTS FOR SYSTEMATIC BIASES 

The validity of any cosmic shear result depends sensitively upon 
the treatment of systematic errors and the control of observational 
biases. Almost all systematic effects, whether they be due to a back- 
ground gradient, astrometric distortions within the telescope or im- 
perfectly corrected PSF anisotropy, act to increase the observed 
correlations between galaxy shapes. All the effects can mimic cos- 
mic shear and the most important task incumbent upon any weak 
lensing survey is to prove that its systematics are controlled to a 
negligible level. As already described in [J2] the astrometric distor- 
tions in WHT have been corrected for, and the basic data reduction 
was performed sufficiently carefully to eliminate most biases. In 
this section, we discuss further tests for other sources of residual 
systematics. 



4.1 E-B decomposition 

The correlation functions can be r ecast in terms of E (gr adient) and 
B (cu rl) modes of the shear field l Crittenden et al. 2000; Pen et al. 
2002). Gravitational lensing is expected to produce only E modes, 
except for a very lo w level of B modes du e to lens-lens coupling 
along a line of sight ISchneider et all20 021. Systematic effects are 
likely to affect both E- and B-modes equally. The presence of any 
non-zero B mode would therefore be a useful indication of con- 
tamination from other sources. 

E- and B-modes correspond to patterns within an extended 
region on the sky. The separation cannot be performed locally, but 
requires the shear correlation functions to be integrated over a wide 




8 [arcmin] 



Figure 9. E-B decomposition of the shear field observed in our WHT sur- 
vey. The points show the measured E (tangential) modes of the shear field. 
As in figure|3] the theoretical curves are predictions for a ACDM cosmol- 
ogy, with Q m = 0.3, Qa = 0.7, T = 0.21 and erg ranging from 0.7 
(bottom) to 1.2 (top). The dashed line line shows our measured B-mode 
(curl) signal, and the shaded region shows its 1<t error bar. In the absence 
of systematics, the B-mode should be consistent with zero. 



range of angular scales. In practice, we cannot perform these in- 
tegrals exactly because our correlation function data extends only 
between ~ 2' and 16'. In other words, a shear field measured within 
a finite aperture can not be uniquely split into distinct E- and B- 
mode components: some correction will always be necessary. 

We find the frequently-used decomposition into aperture mass 
M ap (6>) and M±(6) modes jBartelmann & Schneideife OOO) to be 
particularly unstable with our data. The correlation functions need 
to be integrated between and 28, with a lot of weight placed 
upon their values at small angular scales. Since these are chang- 
ing rapidly, the end result becomes even more sensitive to the bin 
spacing. Furthermore, our measured correlation functions are least 
reliable at small separations. This is likely to be the case in any 
real data because of small-scale effects like overlapping galaxy 
isophotes. Since data from these small scales need to be included in 
all subsequent integrals, any bias there would adversely affect the 
aperture mass on all scales. 

We therefore prefer another method for E/B decomposition. 
Following Crittenden et al. 1 2000), we calculate 



c E (e) = c 1 (e) + 2 



which contains only the i5-mode signal and 

36> 2 \ Ct(#) -C 2 (j?) 



C B (0) = C 2 (6) - 2 



2 



rW , (20) 



which contains only the B-mode signal. These can be calculated 
using data exclusively from scales larger than 8. A correction will 
still need to be made for absence of data on scales larger than 16', 
but this is in a regime where the expected signal (and the neces- 
sary correction) is small. As can be seen from the above equation, 
a function of 9 (not only a constant of integrati on) must be added 
to C E (0) and subtracted from C B {0) (c.f. lPen et alj2002t . We cal- 
culate this function by using theoretical predictions for the best- 
fit cosmological model (as determined in 93.21 to extrapolate our 
data to infinity. The size of this correction is approximately one 
third of the size of the measured iS-mode signal. The correction is 
2.4xl0" 5 at 5' and 1.5 xl0~ 5 at 10'. 

An E/B decomposition of our data is shown in figure [5] On 
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Figure 10. Con'elation functions of stellar ellipticities before (dashed line) 
and after correction (solid points). The shaded region shows the size of la 
errors on the pre-correction correlation functions. Application of the KSB 
method has successfully reduced these by about four orders of magnitude. 



scales 1.7' < 8 < 15', we find a B-mode signal consistent with 
zero, confirming the absence of systematics on these scales. Note 
that because of the extra uncertainty introduced by this additive 
function, we have not used the derived E-mode signal to fit cosmo- 
logical parameters. We instead direcly fit the measurements of C\ 
and C2 on those scales deemed free of systematic errors only (see 

G3. 



4.2 PSF correction 

The WHT PSF over long exposures can be quite anisotropic, with 
a mean stellar ellipticity of 0.051 ± 0.28, where the error quoted is 
the rms stellar ellipticity within one field, averaged over all fields. 
Application of KSB reduces this to 0.0056 ± 0.0012. Figure [Tol 
shows the full correlation functions 



Cf 



( e* e* ) 



(21) 



where i = {1,2}. Application of the KSB method has success- 
fully reduced these by four orders of magnitude. Note that these 
are correlations of ellipticity rather than shear, and therefore ex- 
hibit a different overall normalisation that is connected to the shear 
susceptibility factor. The elimination of stellar anisotropy is also a 
slightly different problem than the challenging correction of galaxy 
shapes, where the interpolation of the PSF is critical, and the size 
of the weight function r g may be different. 

The correction of galaxy shapes for the PSF anisotropy can be 
tested using cross-correlation functions between corrected shears 
7i from galaxies and stellar ellipticity e* before correction, 



( 7> e| ) 2 



(22) 



where i = {1, 2}. These cross-correlation functions are shown as 
dashed line in the top two panels of figure|5| and is consistent with 
zero on all scales. Note that PSF correction residuals would also 
have appeared as B-modes in figure|5| which are in fact consistent 
with zero. 

Uncertainties still remain in the KSB method about the 
overa ll calibration of shear estimators after PSF c orrection (e.g. 
iBernstein & Jarvisl 2002tlVan Waerbeke et alj|2002l). Detaile d im- 



age simulations by B^cone^ lT|2Tjo TrOTlFAeri^rS] 1200 ll) have 
been used to study these issues. iBacon etalJfeOOll) found that, with 




y [arcmin] 

Figure 11. Galaxy shears as a function of position on the WHT CCD, av- 
eraged over all the fields used in our survey. That this is consistent with 
zero everywhere demonstrates an absence of systematics concerning CCD 
readout, telescope vibration and tracking. 




Distance from optical axis [arcmin] 

Figure 12. The radial and "curl" components of galaxy shears as a function 
of distance from the field centre/optical axis, averaged over all the fields 
used in our survey. That this is conststent with zero on all scales demon- 
strates an absence of systematics concerning alt-az tracking, and the suc- 
cessful correction of astrometric distortions within the telescope, whose ex- 
pected behaviour is purely radial (see J2J- 



our survey and telescope parameters, our implementation of KSB 
requires a constant calibration factor of (0.85 ± 0.04) -1 to be ap- 
plied to all calculated shear estimators. 



4.3 Shear as a function of CCD position 

Problems with read noise or charge transfer efficiency on the CCD, 
or telescope flexure and vibration, could cause the measured shear 
to vary as a function of position on the chip, even when averaged 
over many separate fields. Figure fTTI shows plots of shear as a func- 
tion of x and y. Figure fT2l shows plots of shear as a function of r. 
Both are consistent with zero, as desired. 

The mean components of shear across the entire CCD are 
(71) = (1.1 ± 7.1) x 10~ 4 and {72) = (15.6 ± 7.0) x 10" 4 . 
The rms shear within the survey is <r 71 = 0.293 and <r 72 = 0.292, 
or <T| 7 | = (|7| 2 ) 1//2 = 0.413. (Note that our shear measure- 
ment pipeline includes a catalogue cut at |e < 0.5). The main, 
and irreducible, component of this dispersion comes from the in- 
trinsic ellipticities of source galaxies. From other w ork performed 
with h igh resolution and high S/N space- based data iRhodes et alJ 
|2004 iRefreeier Rhodes. & Grothl 12002) and simulated images 
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r g [arcsec] Shear estimator from WHT pipeline [2, 

Figure 13. Left panels: The shear susceptibility factor P 1 for all galaxies, as a function of their size r g , with different options for the signal-to-noise cut in the 
catalogue. The greyscale shows the number density of galaxies throughout the parameter space, calculated by the COMBO-17 pipeline. The WHT pipeline 
produces results almost identical to the top-left panel. In practice the fit is performed on galaxies for one image at a time, and is therefore more noisy; the fit to 
all galaxies here is shown merely to guide the eye. The fit in the top panel is reproduced as a dashed line in the bottom panel. When the COMBO-17 pipeline 
uses a signal to noise cut of v > 5, the population ensemble's average shear susceptibility is raised. Middle panels: A comparison of the shear estimators 
derived via the WHT and the COMBO-17 shear measurement pipelines. Only subsequently matched galaxies are included in this plot, and the WHT pipeline 
always uses a signal to noise cut to the catalogue of v > 15, so all galaxies shown are brighter than that. Consistent results are obtained in the top panel, 
where the COMBO-17 pipeline also uses a cut of v > 15. The unbiased, 1% dispersion of individual values is merely due to the different interpolation of the 
PSF across individual fields. When the COMBO-17 pipeline uses a signal to noise cut of v > 5, the shear estimators of all galaxies are lowered, including 
these bright ones. Right panels: Constraints upon cosmological parameters Q m and erg, from a maximum-likelihood analysis of shears measured with the 
COMBO-17 pipeline, using different options for the signal-to-noise cut. The 68.3%, 95.4% and 99.7% confidence contours are shown. 



lMassevetaI] E004). we estimate a fundamental lower limit for 
<7| 7 | around 0.30. 



5 AN INDEPENDENT ANALYSIS 

Our measurement of as is at the relatively high end of the dis- 
tribution of published cosmic shear results. Compared to equiv- 
alently deep surveys, it is m ost different f rom t hat published 
by the COMBO-17 survey by Bro wn et al] <2003l) . who obtain 
cr 8 (fW0.3) - 49 = 0.721°:m- We have therefore re-analysed our 
WHT data using several variations of the COMBO-17 pipeline to 
determine the extent to which this disagreement might arise from 
technical differences in the pipelines. 

There are four main differences between the COMBO-17 and 
the WHT pipelines. In the COMBO-17 pipeline: 

• the PSF was interpolated across every field via third-order 
polynomials in x and y, so the correction may be different in idi- 
vidual cases. Data was excluded on all scales for any fields with 
problematic PSF correction, 



• a less stringent signal-to-noise cut was applied. All galaxies 
with v > 5 were included, rasising the overall number density 
to 19.8. Every galaxy was still given the same weight when their 
shears were combined, 

• galaxy pairs were not included if the galaxies lay on different 
CCDs. This will lower the signal-to-noise on large scales but may 
remove some bias if the chips are imperfectly aligned on the focal 
plane, 

• no correction was made for astrometric distortions in the tele- 
scope. This will spuriously increase the signal on large scales, 
where galaxies are a long way from the telescope axis, 

• independently-developed code was used to calculate theoreti- 
cal models and to fit cosmological parameters to the data. 

5.1 Shear measurement 

The first two changes are relevant in the creation of a shear cat- 
alogue. Differences in the PSF interpolation introduce a ~1% 
rms dispersion between shear values measured by the WHT or 
COMBO-17 pipelines (see the top-middle panel in figure Il3l This 
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small difference can be explained by the uncertain nature of any in- 
terpolation between sparsely sampled data points. The lower-order 
fit in the WHT method typically produces smoother, and therefore 
perhaps more reasonable behaviour in regions of the images that are 
devoid of stars. However, neither method appears to bias the shear 
estimators or resulting constraints on cosmological parameters. 

Of greater importance are the ^>15ori^>5 catalogue cuts 
in signal-to-noise. Very faint galaxies were discarded in our original 
analys is of WHT images, and the simulated images of Bac on etail 
(2001), because of residual correlation found between the galaxy 
and stellar ellipticities after correction. Such correlations were not 
observed in the COMBO- 17 data, so the cut was lowered. A possi- 
ble explanation is the absence of CCD fringing in the COMBO- 17 
data. Uncorrected fringing in WHT images would alter the shapes 
of those galaxies at a flux level comparable to the fringing via an 
additive process rather than a convolution or multiplication. This 
would not have been corrected, so the galaxies would have sim- 
ply been discarded. Brighter galaxies that appear in both the WHT- 
pipeline and COMBO- 17 pipeline have identical estimates of P 1 , 
since both are calculated by HFINDPEAKS. However, these values 
are noisy. Both pipelines regard galaxies as a population ensemble, 
for which it is possible to average over noise by fitting a global 
shear susceptibility factor. 

Theoretically, the susceptibility factor of a galaxy shear is ex- 
pected to vary as a function of its radial profile, size and ellipticity; 
but not as a function of its flux. We therefore fit P 1 — P 7 (r 9 ), as 
shown in the left-hand panels of figure [T3l However, it seems that 
real faint galaxies do indeed have higher shear susceptibilities than 
similarly-sized bright galaxies. As shown in the middle panels of 
figure H3l this raises P 7 (r 9 ), and thus lowers shear esimators by 
a factor of 0.851° 'g 2 . The skewed distribution reflects the overall 
size distribtuion of galaxies, and arises because this process affects 
small galaxies more than large ones. The reason for this variation 
of P 1 (y) is not yet clear: it may be that the morphology distribu- 
tion of faint galaxies really does contain intrinsically higher shear 
susceptibilities. Both approaches would be valid in this case, since 
the fit really has found a suitable average value for the population 
ensemble of galaxies. 

A more worrying alternative would be that the noise in faint 
galaxies biases P 1 to hig her values, partly because the ma in P Bh 
component of P 1 (see lKaiser. Squires & Broadhurstll 19951) is de- 
fined to be strictly positive. In this case, it would instead be prefer- 
able to fit only the bright galaxies, then apply that susceptibility 
to faint galaxies of the same size. Either way, the results from the 
WHT pipeline would be less affected, as the v > 15 cut has been 
calibrated upon simulated images containing a known signal by 
iBaconetalJfcOOll) . 

A full investigation into this technical issue is beyond the 
scope of this paper, and it will need further investigation in the fu- 
ture. Other comparative studies, e.g. the Shear TEsting Program 
(STEP, Heymans et al. in preparation), or the Edinburgh and Bonn 
pipelines (Heymans et al. in preparation), find that the detailed way 
in which P 1 is, or is not, adequately fit may indeed explain a large 
part of the variation between results from different implementations 
of KSB. Here, we shall propagate the analysis of shear catalogues 
from the COMBO- 17 pipeline with both v > 15 and v > 5 cuts. 



5.2 Cosmological parameter constraints 

We find similar instabilities in the production of shear-shear corre- 
lation functions with the COMBO-17 pipeline as we found in the 
WHT pipeline; results are sensitive to binning at the ~ ±7% level. 



Measurements with the v > 5 as shown as blue circles in figure|3| 
rescaled by a factor (0.8/0.9) 165 to allow comparison despite the 
increased source redshift distribtion. Excluding galaxy pairs of dif- 
ferent CCDs has increased the size of error bars by up to 50% at 
large 9, where the number of available galaxy pairs is significantly 
lower, but the central values move only within their la error bars. 
Not correcting for the telescope's known astrometric distortion has 
moved the bin at largest scales by a significant amount, but this is 
not used for parameter fitting anyway. Excluding data on all scales 
from fields with any residual star-galaxy correlation also has a min- 
imal effect. Each of these differences change the derived constraint 
on ag by less than 1%. 

The main difference in the two analyses is that introduced 
by the different cuts in signal-to-noise. The inclusion of fainter 
galaxies increases the median magnitude of the v > 5 popula- 
tion by 0.84 ± 0.4 mag nitudes, leading to a new median redshift 
of z s = 0.90 ± 0.07 iCohen etalJ <2000l) ). The expected cos- 
mic shear signal thus increases (see equation 1 14i - an effect that 
is incorporated during the calculation of theoretical models for the 
fitting of cosmological parameters. Reassuringly, the shear signal 
measured from faint galaxies is also higher. The final constraints 
for the COMBO-17 pipeline, using a cut of v > 15 indicate: 



an 



03 



1.01 ±0.13 , 



(23) 



including only statistical errors and, for a cut of v > 5 with the 
higher assumed source redshift distribution 



/n m \°- 52 

CT8 V"oTJ =0.98 ±0.10. 



(24) 



The consistency of these fits, and their agreement with the 
results from the WHT pipeline is reassuring. They tend to justify 
the adoption of a global shear susceptibility factor for an ensemble 
population of galaxies, and verify the validity of the methods for 
controlling many potential systematic effects in both pipelines. 

Our survey probes a wide range of angular scales, and thus 
excludes small values of f2 m < 0.25 (at 68%CL). The prefer- 
enc e for large f2 m com es from the fit to the transfer fun ction 
by jBardeen et aljfl986l) . If we adopt the fit by lEisenstein & HuJ 
1 1997), and assume the Hubble parameter h — 0.7, we ef- 
fectively apply a different prior on Q m . In this case, our data 
excludes Q. m > 0.55 (at 68%CL), similar to th e result ob- 
tained by IVan Waerbeke. Mellier & Hoekstri] (2005). Using the 
lEisenstein & Hul 1 19971) ' transfer function also lowers our best-fit 
value for as around S7 m = 0.3 by ~2%. 



6 CONCLUSIONS 

We have measured the weak lensing shear-shear correlation func- 
tions in four square degrees of deep ii-band imaging data from the 
William Herschel Telescope. Our measurements constrain the am- 
plitude of the mass power spectrum, cr8(fim/0.3) ' 52 = 1.02 ± 
0.15, including all contributions to the total 68%CL error budget: 
statistical noise, sample variance, an additional estimated error due 
to binning instabilities from non-Gaussian outliers, covariance be- 
tween different angular scales, systematic measurement and detec- 
tion biases, source redshift uncertainty, and marginalisation with 
priors over other parameters. We have examined our data for con- 
tamination by systematic effects using a variety of tests including 
an E-B decomposition. These demonstrate a well-understood and 
modest contribution to our uncertainties from systematic errors. 
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Using the pipeline developed for earlier WHT studies, we 
find our measurement of the normalization of the dark matter 
power spectrum lies at the relatively high end of the distribution 
of published values. However, it is still consiste nt with those from 
equi valently deep surv eys by R efregier Rhodes. & Groth I <2002l) 
and iRhodes et all2004l) . Our results are also consistent at the la 
level with CMB results from the Wil kinson Microwave Anisotropy 
Probe (WMAP) ISpereel et all2003l) . 

The wide distribution of a$ constraints from recent cosmic 
shear surveys understandably casts some aspersion upon their pre- 
cision. For example, it might be argued that the dispersion largely 
arises from unknown or poorly-understood systematic effects. For 
the first time we analyze our dataset with an independent pipeline 
- that developed for the COMBO- 17 data <Brown et all J2003)). 
This provides a valuable check on the extent to which dispersion in 
the published cosmological results, as well as our apparently high 
normalization, might arise from different techniques, both in con- 
structing shear catalogues and in analysing the shear-shear correla- 
tion functions. Reassuringly, we find remarkable concordance be- 
tween the two pipelines for the same dataset. However, the various 
differences we have explored, for example in the selection thresh- 
olds, seem to be insufficient to reconcile the spread in observed as 
values, suggesting that significant differences remain at some level 
in the actual data themselves. 

Uncertainties in the redshift distribution of source galaxies in 
deep data clearly contribute. It is difficult to determine the precise 
redshift distribution of galaxies after excluding those smaller than a 
fixed apparent size. We have been conservative in this analysis and, 
as seen in equation 1151 . source redshift uncertainty is already a ma- 
jor component of our total error budget. The resolution of such is- 
sues will require extensive spectroscopic follow-up and more com- 
plete image simulations. Such advances are essential if the poten- 
tial of the next generation of cosmic shear surveys is to be fully 
realised. 

Finally, we note that most recent cosmic shear results re- 
main discrepant at the 3a level with measurements derived 
from the abundance of X-ray selected cluster samples based 
on an observational rather t han theoretica l mass-temperature re- 
lation iBorgani et alJ 1200 it ISeliak j [20011 iReiprich & Bohringed 
| 2002t IViana. Nichol & Liddldl2002l> . These suggest a s « 0.75. 
Amara & Refregied J2004h concluded that even extreme non- 
Gaussianity in the mass distribution would be insufficient to ex- 
plain this discrepancy, because the two techniques probe simi- 
lar mass scales. Further studies are therefore needed in both the 
cluster method, to understand the difference between the observed 
mass-temperature relation and that found in numerical simulations; 
and in the weak lensing method, to construct more reliable and 
better calibrated shear measurement methods. Such consistency 
checks will represent a crucial verification of the standard ACDM 
paradigm, so resolving this issue is of paramount importance. 
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